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ABSTRACT 

A new improved model for cosmic rays - middle atmosphere 
interaction is developed. The ionization q(h> -profile dependence 
on penetrating high energy particles composition (protons, 
alpha-particles and heavier nuclei) and energy spectra (solar 
activity modulation included) are investigated. A computer 
program, realizing the Gaussian algorithm for solving of 
multidimensional integrals is created. The correspond ing electron 
density profiles N(h) at solar minimum and maximum are obtained, 
too. 


INTRODUCTION 

The galactic cosmic rays (GCR) flux is an important 
factor for the middle and low atmosphere ionization under day and 
night conditions. GCR form the independent CR— or C— layer at 
height from 50 to 80 km in the ionosphere, thus forming its base 
(VELLINOV et al., 1974). Here - at 50-60 km height a region 
char acter i zed with anomalously high field values up to — 4f— 6 V/m 
(APSEN et al., 1988; KQCHEEV et al. , 1976, MAYNARD et al«, 1981) 
is registrated via rocket measurements of the atmospheric electric 
field E^-prof iles. Thus, GCR ionization proves important for the 
clarification of a number of problems concerning the middle 
atmosphere electrodynamics and middle atmosphere — ionosphere 
interaction. For that purpose, a new model for the GCR - middle 
atmosphere interaction will be developed that will contribute to 
the further precisian and generalization of previous results 
(VELLINOV et al., 1974). 


GCR IONIZATION IN THE MIDDLE ATMOSPHERE 

The electron production rate (cmTs. > at height h (km) for 
the particles type i is 

f ' ^:0 1 

where Q. =35 eV Is the energy necessary for one electron-ion pair 

formation, (dE/dh) - the particles ionization losses, D(E) - their 
differential spectra,^- the azimuth, 0- the angle towards the 
vertical. As GCR penetrate isotropically from the upper 
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where E- is 1:he energy (GeV/nucl) correspond ing to the geomagnetic 
threshold R(GV>, A0 takes into account that at a given altitude 
the particles can penetrate from spatial angle j0,90 +A0) that is 
greater than the upper hemisphere angle - (0,90 # ) . For AS there 


is the following equation (VELLINOV, 1960): 

A 8=90* -arc cos(2R # h+hV^R 4 +h) <3) 

where R =6371 km is the Earth radius. In the interval 30-100 km 
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0 changes from 6 to 10*. 

For dE/dh we are going to use the dependency 

(DOBROTIN, 1964; VELLINOV et al. f 1974, p.202) 

j ^ MtV. E“ ]0 -f (l.fi GeV region 1 (4) 

f MzV. £*' . £* O.G-^/0~ 3 &cV region 2 (5) 

where j> is the atmosphere density in g.cfiT 3 . In fact the 
particle energy fall law in region 1 (0,6-10 GeV) is 

E(h)=E - 2.10’ 3 hZ X Ch <0 , h) (6) 

r* r -i 

where h= J o(h)dh is the depth (g/cm ) , Z is the particle charge, 
and Ch(0,h) Chapman function taking into account the spherical 
form of the atmosphere. During particles penetration in the 
middle atmosphere their energy decreases while their ionization 
lasses in the region 1 remain constant; in the region 2 they rise 
drastically according to (5): That effect was not taken into 

consideration in the previous investigations. After deducting 
(4,5) in (2) the total equation for the electron production rate 
is received: 


q(h)=q 1 (h>+q jt (h) 


(7) 


where q lJ from the region 1 practically coincides with the 
previous calculations ( VELLINOV , 1968) - The ionization in the 
region 2 will be added to q. 
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In fact, the function 
integral spectra in the regions 1 
atomic weight. The angle 8, 


f , (8) 
and 2 


equation (10) , when the proportion TY* <8 


normalizes the particles 
A is the particles 
determined with the help of 
h ) =K ( h ) is fulfil led. 




COMPUTER REALIZATION OF THE MODEL AND ANALYSIS OF THE 

RESULTS 

The theoretical model thus presented is realized with the 
help of a computer. The Gaussian method is applied for the 
solution of multidimensional integrals (PRESS et al., 1987) . It is 
character i zed by high precision for the smooth integrant 
functions, that are well approximated with polinoms. That method 
is machine run— time efficient — it uses only 10 function values in 
the integral interval. 
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Having in mind, that Z/A=l for the protons and about 0,5 
for the heavier nuclei, as well as the actual GCR composition 
(DORMAN, 1963) q.(h) were calculated for all groups of nuclei* 
protons (Z=l), alpha-part icl es <Z=2), light <Z=3-r5), middle 
<Z=6 t 9>, heavy (Z = 10-r21) and very heavy (Z>22>. In consequence 
the ionizations of the separate groups were summed up and the 
total ionization q(h) was received. 

The calculations were conducted for two GCR spectra - in 
maximum and minimum solar activity, i.e. the 11-year GCR 
modulatiion from the solar wind is considered (DORMAN , 1963) . That 
modulation is expressed in the reverse solar activity changes and 
the GCR flux. A part of the results from the program model 
realization are presented on the Table: q(h) in the upper (B0 km) 

and lower (50 km) part of the C— layer. The results from the 
previous calculations are given in brackets (VELLINOV et al. , 
1974). Hence, the contribution of to the total q can be 

considerable. From the Table it is seen also that at geomagnetic 
latitudes over 55° q(h) is growing during the solar minimum 
period two times, compared to the solar maximum. At middle 
latitudes the 11-year variation is still considerable, while at 
low and equatorial latitudes the variation is weak. 


C-LAYER ELECTRON DENSITY DURING MAXIMUM AMD MINIMUM SOLAR 


ACTIVITY 

equation 


The equilibrium electron density (cm 


-3 , 


JfW- 




1 V* 




is obtained from 


where ^ =N~/N is the negative ions-electrons ratio and is the 
effective recomb inat ion coef icient. Having in mind the values 
from MITRA (1986) the C-layer electron density was calculated for 
typical day conditions for middle (41*) and high (55 # ) geomagnetic 
latitudes during solar minimum and maximum. The results are given 
on the Figure and are mean valued in relation to season dependency 
due to the unsufficient representativeness of the seasonal 
variations- But, it can be said definitely, that in the C— layer 
(and below it) maximum in winter N is 30—40% greater than in 
summer. Over 80 km the season variation changes its sign. 


CONCLUS I ON 

The present work consideres new aspects of the C-layer 
which together with the total GCR ionization in the middle and low 
atmosphere is of major importance for the global 
atmospher ic-electr ic circuit in the system of the 
solar— terrestr ial relationships and biorelationships. 

The results obtained so far are not final and are an 
object of further specification. For example under significant 0 
the particles pass greater quantity of substance connected with 
the drastic growth of the Ch function. In that case an object of 
consideration are not only the electromagnetic interactions of GCR 
(4,5) but the nuclear as well- And those secondary cosmic rays 
will probably increase still more the obtained herewith electron 
production rates. 

Under high solar activity the importance of the solar 
cosmic rays in ionization will grow, and hence — in the 
conductivity, currents, the electric fields and energetic 
processes in the middle atmosphere (VELLINOV and MATEEV, 1989). 
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C-Layer Electron Production Rate q.cm s Table 

HIBH LATTITUDES - ABOVE 55*, R=1.5 GV 


h , km 

SOLAR MINIMUM 

SOLAR MAXIMUM 


SUMMER WINTER 

SUMMER WINTER 

80 

7.5 10" J 3,7 10" J 
C5 10 2.5 10* 3 

4 10“ J 2 10 -3 

2.6 1C' 3 1.3 10' 3 3 

2.5 10"' 1.3 10*' 

50 

4,8 10" 2.8 10“' 


C2.5 10*' 1.5 10"' 

1.3 10*' 7 10"' 3 


v . = >= 

MIDDLE LATITUDES 41* R=5 GV 


ff 80 

3.6 

10' J 

1.8 10' -1 

2.7. 10" 3 

1 .5. 10‘* 


C2.4 

10“ J 

1.2 10' 3 

1 .8. 10“ s 

10" 3 3 

[50 

2.3 

10“' 

1.3 10”' 

1.7. 10“' 

10“' 

1 

1 

Cl .2 

10“' 

7 10 

9. 10 

5.5. 10“* 3 


h , km 



latitudes during maximum and minimum solar activity. 
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